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ABSTRACT 

One explanation for the enhanced ratio of volatiles to hydrogen in Jupiter’s atmosphere compared 
to a a gas of solar composition is that the planet accreted volatile-bearing clathrates during its for¬ 
mation. Models, however, suggest that S would be over abundant if clathrates were the primary 
carrier of Jupiter’s volatiles. This led to the suggestion that S was depleted in the outer nebula 
due to the formation troilite (FeS). Here, this depletion is quantitatively explored by modeling the 
coupled dynamical and chemical evolution of Fe grains in the solar nebula. It is found that disks 
that undergo rapid radial expansion from an initially compact state may allow sufficient production 
of FeS and carry H 2 S-depleted gas outward where ices would form, providing the conditions needed 
for S-depleted clathrates to form. However, this expansion would also carry FeS grains to this region, 
which could also be incorporated into planetesimals. Thus for clathrates to be a viable source of 
volatiles, models must account for the presence of both H 2 S in FeS in the outer solar nebula. 

Subject headings: planets and satellites: formation - planets and satellites: composition - planets and 
satellites: individual (Jupiter) - protoplanetary disks 


1. INTRODUCTION 

In the core accretion model, Jupiter’s atmosphere is 
largely derived from solar nebula gas that was gravita¬ 
tionally captured by a ~10 Earth mass solid core, with 
additional contribution from planetesimals that were ac¬ 
creted during the plan ets formation ( [Pollack et al.| 1996 
Hubickyj et al.||2005|). As the nebular gas would have 
had roughly a solar composition, any deviations in the 
composition of the planet’s atmosphere from that of the 
Sun is thought to reflect the contribution of the late ac¬ 
creted planetesimals. The Galileo probe, during its de¬ 
cent, found that the the ratios Ar/H, Kr/H, Xe/H, C/H, 
N/H, and S/H we re 2-4 times greater than a gas of solar 
composition gas (Owen et al. 1999). The high abun¬ 
dances of these volatile elements, particularly the noble 
gases, is puzzling as they do not react with solids and do 
not freeze out as pure species, and thus are not expected 
in planetesimals until very low temperatures (<~10 K) 
are reached. 

Two mechanisms that attribute the enhanced abun¬ 
dances of these elements to accretion of cold, icy plan¬ 
etesimals have been proposed. The first was that these 
elements were trapped in amo rphous water i ce at very 
low temperatures (<30 K) (Owen et al. 1999). At these 
low temperatures, experiments have shown that gaseous 
molecules and atoms that are adsorbed on a solid sur- 
face become trapped as water freezes out on that solid 


dLunine fe S tevenson|1985[|Gautier et al.|2001|). [Gautier] 
et al.| (|200ip calculated the abundances of volatiles that 
would be trapped as clathrates formed in a cooling so¬ 
lar nebula and found that the predicted abundances of 
noble gases, C, and N would match the abundances in¬ 
ferred from the Galileo probe measurements. However, 
they found that S would be present at a greater abun¬ 
dance than was observed. 

To explain the discrepa ncy in the predict ed and ob¬ 


served S abundances, Gautier et al. (2001) suggested 


that H 2 S, the primary gaseous carrier of S in the planet 
formation region was depleted due to the formation of 
troilite via the reaction: Ee( s )+Ff 2 5'( s ) —> FeS^-\-H 2 [ g ). 
Troilite is commonly found in meteoritic materials, and 
the reaction occurs on timescales that a re sh ort compared 


et al. 1996 


Tachibana & 


to the lifetime of the sola r nebula (Fegley|l988 Lauretta 

Tsuchiyama 1998]). |Gautier 
ly ^35% of H 2 S would 


(Bar-Nun et al.||1985[ |Notesco et al.||2003[ |Yokochi et al.| U veiy unsu ggortable. 
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et al. (2001 argued that only ~35% of H 2 S would need to 
be removed via this reaction to match the Galileo probe 
measurements, and suggested that as FeS formed in the 
terrestrial planet region of the solar nebula, H 2 S depleted 
gas diffused out to the Jupiter formation region, allowing 
for lower amounts of S to be locked up in clathrates than 
would otherwise occur. While some use such a depletion 
for a starting point i n giant planet format ion models (e.g. 
Mousis et al.||2014|), |Atreya et al.| (|2003|) challenged this 
idea, calling the removal of this much H 2 S “quantita- 


2(JI2 H As very low temperatures are required to pro- 
duce amorphous ice and to allow for efficient trapping, 
this mechanism requires that Jupiter formed from or ac¬ 
creted planetesimals from much further away from the 
Sun than its current location, or that temperatures in 
the Jupiter-formation zone were much lower than pre¬ 
dicted by current models for the solar nebula. 

The second mechanism involves forming clathrate hy¬ 
drates, crystalline ice with guest molecules trapped inside 


Lodders (2004) offered an alternative explanation for 
the volatile abundance. Citing the low oxygen abundance 


measured by the probe, Lodders 
formed inside the snow line - 


(2004) suggested Jupiter 
1 volatiles were then ac¬ 


creted directly from the solar nebula gas, with additional 
C and N delivered in the form of organics. In this model 
sulfur could trace the contribution of solids and gas ac¬ 
creted by the planet: if S was distributed purely in the 
gas, then its relative abundance should be the same as all 
other gaseous species; if S was present as a solid (FeS) 
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then its abundance would reflect the amount of solids 


accreted by Jupiter as it formed. Lodders (2004) noted 
that the S-to-noble gas ratio was solar, and assumed that 
it was delivered as FeS, and thus concluded that Jupiter 
formed from a solar mix of solids and nebular gas. 

The behavior and speciation of S in the solar nebula 
and its distribution between the gas and solids is thus 
critical to understandin g how volatil es were delivered to 
a young Jupiter. Pasek et al. (2005) explored this issue 
by modeling the evolving concentration of H 2 S within 
a turbulent, diffusive protoplanetary disk. H 2 S was as¬ 
sumed to be removed from the gas due to FeS formation 
where temperatures fell in the range of 500-690 K. The 
upper bound of 690 K was chosen as this was where FeS 
became stable in the nebula, while the lower bound of 
500 K was set because reaction kinetics were too slow 
to produce much FeS. Any F[ 2 S that diffused into the 
region of the nebula that fell in this temperature range 
was assumed to immediately react with Fe and the sulfur 
permanently locked up there as FeS. As a result, the FI 2 S 
abundance decreased throughout the nebula over time, 
reaching <65% of its i nitia l value inside ~6 AU in 10 5 
years, as predicted by |Gautier et al.| (f200Tj) . 

While providing important insights into the F[ 2 S evolu¬ 
tion in the solar nebula, th e treatment o f FeS f ormation 
was somewhat simplified in Pasek et al. (2005). For ex¬ 
ample, chemical kinetics were calculated assuming the 
pressure and temperature conditions at the disk mid¬ 
plane and a solar abundance of H 2 S. Due to diffusion, 
grains would have been mixed throughout the vertical ex¬ 
tent of the disk, where temperatures and p ress u res wer e 
much lower, and reaction timescales longer (Ciesla|2010). 
Further, the lower abundances of H 2 S as the disk evolved 
implies longer reaction times as the rate of FeS formation 
is directly proportional to the partial pressure of H 2 S. Fi¬ 
nally, particles would have migrated in and out of the FeS 
formation region on timescales less than or comparable 
to the chemical timescales, suggesting that complete re¬ 
action between H 2 S and Fe would not necessarily occur. 

Here a new model that tracks the motions of dust 
grains in a dynamically ev olving proto planetary disk is 
used to expand the study of [Pasek et al.|(2005]) and quan¬ 
tify in detail the extent to which Fe and S would react 
prior to Jupiter’s formation, with the focus on the level 
to which H 2 S becomes depleted as a result of troilite for¬ 
mation and the implications for the S abundance where 
clathrates could form. The modeling methods are de¬ 
scribed in Section 2. In Section 3, the results for two ex¬ 
amples of protoplanetary disk evolutionary histories are 
presented. These results and the implications for the de¬ 
termining factors for planetary compositions in our Solar 
System are discussed at the conclusion of the paper. 


2 . 1 . 


2. MODEL 

Particle and Disk Evolution 


P article tracking methods 


& Sandford 


2012 Charnoz et 


(Ciesla] 

aT 201 


2010 2011 Ciesla 


I) are employed to 


determine the exact physical environments the Fe/FeS 
grains were exposed to within the evolving solar nebula. 
The evolution of the protoplanetary disk is modeled using 
the classical a-viscosity formalism, where stresses within 
the differentially rot ating gas drive mass a nd angular mo¬ 
mentum transport (Shakura & Sunyaev 1973 Lynden- 


Bell & Pringle 1974 Hartmann et al. 1998| ), which al¬ 
lows us to consider the evolution of the disk over ~T0 6 
years of evolution. The viscosity in the disk is given by 
v=ac s H 1 where c s is the local speed of sound, H is the 
local disk scale-height, and a is a parameter describing 
the strength of the viscosity. The motions of the grains in 
the disk were calculated accounting for vertical settling 
by gravity, gas drag migration, advective flows associated 
with disk evolution, and diffusion, where the diffusivity 
was taken to equal to the disk viscosit y {V ~ v). A 
similar model was used by Pasek et al. (2005), though 
only diffusive redistribution of the gas was considered; 
the disk structure was not allowed to evolve during the 
time of interest. 

Iron particles of a given radius, a, were suspended 
across the disk at constant Fe/H 2 mass ratio of 
e=1.67xl0 -3 ( Palme et al.|2014 ). The mass of the parti¬ 
cles was set at m gra in=^cL i PFe where pFe= 7.874 g/cm 3 . 
As tracking the dynamics of each iron grain in the solar 
nebula is computationally prohibitive (a Mdi S k = b.lMQ 
disk would contain ~1.67xl0~ 4 M 0 of Fe, which would 
amount to ~10 4 ° Fe grains 1 pm grains), each model 
grain is treated as a super-grain with its evolution repre¬ 
senting the evolution of a number of other grains equal 
to N rep = eM disk /(m grain n part ), where n part is the num¬ 
ber of super-grains considered in the simulation. The 
super-grains were distributed radially within the disk to 
ensure a solar-ratio of Fe/H 2 throughout. Grains were 
initially placed at the disk midplane, however vertical 
mixing leads to gra ins becoming well mixed vertically 
on short timescales (Ciesla 2010), making the choice of 
initial Zq unimportant for our results. 

2.2. Gas-Grain Chemistry 

Iron grains sulfidize through reaction with H 2 S to pro¬ 
duce troilite, while troilite decomposes following the re¬ 
verse reaction: FeS( s )+H 2 ( g ) —> Fet^+H^S^y Because 
of the low pressures expected in protoplanetary disks, 
one of the primary controls on the rate of gas-grain reac¬ 
tions is the collisio n rate of rea cting gas molecules with 
the grains (|Fegley| |1988|, |2000|). The collision rate of a 
gaseous species, 1 , with a gram is: 


7 i = 2.5 x 10 : 


,31 


Pi 

MiT 


47tot 


( 1 ) 


where P* is the pressure of the species in bars, M, is 
the molecular weight of the gaseous molecule, T is the 
temperature of the gas and solid, and a is the radius 
of the grain. The fraction of such collisions that would 
result in the reaction proceeding would be determined 
by the activation energy of the reaction, E, and given 

E 

by f = e~ ft where R is the gas constant. Experiments 
found that the activation energy for troilite formation 
(the for ward reaction here) falls in the r ange of 42-50 
kJ/mol (|Pasek et al.||2005| La,uretta||2005[). We adopt a 
value of Ej=5 0 kJ/mol, but have found lower values in 
this range do not change the results significantly. The 
reverse reaction, where the S in FeS combines with H 2 to 
form H 2 S gas leaving Fe behind, can be determine d using 
t he eq uilibrium constant for the reaction as in Lauretta 
et al. (1996), giving E r =115 kJ/mol. The rate at which 
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iron is converted into FeS is given by: 


dNpe 

dt 


= -4tt a * 2 3 * (7 H 2 sff/3 - 7 H 2 fr (1 - 0)) (2) 


where Np e is the number of Fe atoms in a grain and /3 
is the fraction of iron in the grain remaining in metallic 
form, while (1 — /3) represented the fraction of iron taken 
up as FeS. Only small, micron-sized grains are considered 
so that the diffusion of S in the grai n is not a rate-l imiting 
step, and thus can be ignored (Pasek et al. 2005). 


2.3. Calculation Method 

Two illustrative cases are considered, each with a 0.1 
M 0 disk with an initial surface density profile given by 

E (r)=E 0 extending from r=0.05 AU to r max , 

with r max =10 and 50 AU at the beginning of the simu¬ 
lations. In each case a=10 -3 . Grain movemen t is calcu ¬ 
lated using the particle tr acking methods in Ciesla fe 


Sandford 


(2012) (see also Ciesla 2010 2011 Charnoz 


et al.||20ll ). In carrying out the calculations, the disk 

evolves for a time step, yr, after which time 

the surface density and temperature are updated every¬ 
where. The movement of the grains are then determined 
using the local disk properties. The chemical evolution 
of these grains is calculated using the conditions at the 
new location. In calculating the chemical evolution of 
each grain, it is assumed that the solids remain in contact 
with the same gas throughout the course of a calculation. 
That is, as troilite forms in a grain, that grain continues 
to react with the H 2 S-depleted gas. Grains that migrate 
within r= 0.5 AU of the Sun are assumed to be lost from 
the disk, carrying whatever S they contain with them. In 
all cases, n part = 2xl0 5 super-particles are considered in 
the model, with each grain representing ~5xl0 34 actual 
iron grains. 

The vertical structure of the disk is calculated assum¬ 
ing a purely viscous model as the focus is on the early 
evolution of the disk when mass and angular momentum 
transport are rapid and thus internal heating would dom¬ 
inate. High t emperatures are expected during this early 
epoch (Boss 1998). At later times, mass transport may 
be limited to the surface layers of the disk and irradia- 
tive heating would control the thermal structure leading 
to temperatures throughout the disk would be so low to 
prevent significant reaction between H 2 S and Fe or de¬ 
composition of FeS. This meaning that iron sulfidization 
was limited to the very early stages of solar nebula evo¬ 
lution. Note that the changes in H 2 abundance due to 
chemical reactions are minimal and do not impact the 
overall surface density of the disk anywhere. 

3. RESULTS 

Figure 1 shows the results of the two cases described 

above, with the H 2 S distribution in the gas, normalized 
to its solar value, shown throughout the inner 10 AU of 
the disk at times of 0.25, 0.5, 0.75 and 1 Myr. The dis¬ 
tribution is found counting the particles at each time in 
a series of radial bins and calculating the amount of H 2 S 
remaining in for each particle, then dividing by the area 
to get the surface density of the gas. This is then nor¬ 
malized to the abundance of H 2 S that would be expected 
for a solar S/H ratio, given the local surface density. 


Valleys in the H 2 S distribution are readily apparent, 
centered in the 0.5-1.5 AU range in each case. These 
valleys correspond to where FeS was readily produced 
through reaction of Fe with H 2 S, as midplane tempera¬ 
tures ranged here from ~400-700 K. Viscous disks begin 
in a very hot state due to the high rates of dissipation 
(heat production) and high optical depths (heat reten¬ 
tion), but cool with time as they lose mass to their cen¬ 
tral stars and expand due to angular momentum conser¬ 
vation, thinning the disk and allowing radiation to escape 
more readily. Thus the region of FeS production migrates 
closer to the star with time. Inside of these valleys tem¬ 
peratures were higher leading to decomposition of FeS 
and the return of H 2 S to the gas, while in the outer disk 
temperatures were too low for much FeS to form and any 
FeS found there, or H 2 S removed, is a result of reactions 
in the inner disk and outward transport through the disk. 
The gradual decrease in H 2 S near the inner edge of the 
disk is an artifact of the calculations and arises due to 
the removal of particles from the simulation inside of 0.5 
AU to represent loss by accretion onto the star. 

While similarit ies exist bet ween the results presented 
here and those of Pasek et al. (2005), in that H 2 S abun¬ 
dances can decrease in the disk over time, key differences 
are apparent. Most important among these is that the 
concentration of H 2 S do e s not decrease indefinitely here 
as it did in Pasek et al. (2005). Instead, the amount of 
H 2 S outside the FeSTormation valley remain relatively 
constant for the r max =50 AU case. While outward dif¬ 
fusion carries H 2 S depleted outward in the disk, the net 
motions of material due to viscous evolution in the disk 
are inward in this region. Thus the H 2 S-depleted gas can 
only diffuse outward so far before being pushed inwards 
by the advective motions associated with disk evolution. 
However, the initially more compact disk (r max =10 AU) 
experiences very rapid outward expansion in this region 
of the disk, meaning that H 2 S-depleted gas not only dif¬ 
fuses outwards, but is also carried outward by the net 
motions of the disk. This leads to a significant overall 
depletion of 2-10 AU region shown here. Efficient de¬ 
pletion occurs only very early in the history of the disk 
however, as once the disk expands the net motions of 
the 2-10 AU region become directed inward. This allows 
materials which originated at greater distances, and thus 
did not become depleted in H 2 S, to be carried inwards. 
As a result, the H 2 S abundance increases with time. 

4. DISCUSSION 

Depletion of gaseous H 2 S in the solar nebula as a re¬ 
sult of troilite formation was possible in the solar nebula. 
The greatest depletion would have occurred in the inner 
nebula, where physical conditions allowed the reaction 
to occur on short timescales. Depletion levels of 65% or 
greater in the o uter s olar nebula as needed in the clathra- 
tion model of Gautier et al. (2001), require that the disk 


begin in a massive, compact configuration. This com¬ 
pact structure would have allowed FeS to form in the 
inner disk with H 2 S-depleted gas being carried outward 
during the radial expansion of the disk that occurred dur¬ 
ing its early viscous evolution. Such a depletion is more 
difficult to achieve in disks which begin with a more ex¬ 
tended structure, as the net flow for the inner disk is 
directed towards the star. 

An initially compact disk would result in significant 
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outward radial transport of mater i als be yond just Fe 


and F eS gra ins or F[ 2 S gas. Ciesla (2010) and Jacquet 
et al. (20111 argued that a similar configuration would 


aTdTn the outward transport and preservation of Cal¬ 
cium, Aluminum-rich Inclusions (CAIs) found in mete¬ 
orites. Indeed, this dynamical evolution could lead to 
significant mixing of inner solar nebula materials into 
the ou ter solar nebula as postulated by Prinn fe Fegley| 
(1989) and be consistent with the high abundance of re¬ 


fractory materi als found in the Stard ust samples from 
Comet Wild 2 (|Brownlee et al. 20061. Thus this evo¬ 
lution may be consistent with other evidence for radial 
mixing and transport in the disk. 

In the calculations presented here, H 2 S was depleted 
only through reaction with Fe to produce FeS, however 
other reactions may occur, incorporating S into solids. 
For example, H 2 S may react with NH 3 gas at lower tem¬ 
peratures to form NH 4 SH solids. The importance of this 
reaction depends on the extent to which N is present 
as NH 3 : equilibrium chemical models predict that NH 3 
should be the dominant form, but sluggish kinetics may 
keep much of N in the form of N 2 (Lewis & Prinn||l980). 
The exact level of conversion will depend on the P-T 
conditions that gas was exposed to during solar nebula 
evolution. This issue will be the focus of future work, 
using the methods outlined here. 

While H 2 S can be depleted in the outer solar nebula, 
this region is not predicted to be depleted in total S. FeS 
grains would be carried outward by the same dynamical 
evolution which carried H 2 S-depleted gas to this regi on. 
Taylor et al. (20041 criticized the model of Gautier et al. 


( 2001 ) on the grounds that removal of the needed sulfur 
in the gas phase would lead to too much FeS in the inner 
solar nebula, which cannot be reconciled with the ob¬ 
served compositions of the terrestrial planets. Here, we 
find that FeS would not remain in the inner solar nebula, 
but would be mixed and made available to outer Solar 
System bodies. As a result, the S content of planetesi- 
mals would be higher than predicted by clathrate models 
alone: while a lower abundance of H 2 S allows the p roper 
S/H r ati o to be achieved in the ices described by Gautier 


et al. (2001), additional S would be found in the rocky 


component of the dust that these planetesimals formed 
from. This suggests that H 2 S must have been depleted to 
an even greater degree if clathrate-bearing plantesimals 
are to properly account for the compositions of Jupiter’s 
atmosphe re. It is unlikely all S is locked up in rocky ma¬ 
terials, as[Boissier et al. (2007) found that H 2 S emissions 
from Comet Hale-Hopp were consistent with the molecule 
being present in comet itself, and not the product of pho¬ 
tochemistry in the coma. Additional measurements of S 
abundances and its speciation in comets will provide im¬ 
portant insights into how this element was distributed 
and the extent to which it may have been locked up as a 
rocky mineral during the early, hot phase of solar nebula 
evolution. 
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